Generalized Slow Roll for Large Power Spectrum Features 
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We develop a variant of the generalized slow roll approach for calculating the curvature power 
spectrum that is well-suited for order unity deviations in power caused by sharp features in the 
inflaton potential. As an example, we show that predictions for a step function potential, which has 
been proposed to explain order unity glitches in the CMB temperature power spectrum at multipoles 
£ = 20 — 40, are accurate at the percent level. Our analysis shows that to good approximation there 
is a single source function that is responsible for observable features and that this function is simply 
related to the local slope and curvature of the inflaton potential. These properties should make the 
generalized slow roll approximation useful for inflation-model independent studies of features, both 
large and small, in the observable power spectra. 



I. INTRODUCTION 

The ordinary slow roll approximation provides a 
modcl-indcpcndcnt technique for computing the initial 
curvature power spectrum for inflationary models where 
the scalar field potential is sufficiently flat and slowly 
varying. Such models lead to curvature power spectra 
that are featureless and nearly scale invariant (e.g. [l[). 

On the other hand, features in the inflaton potential 
produce features in the power spectrum. Glitches in the 
observed temperature power spectrum of the cosmic mi- 
crowave background (CMB) [2[ have led to recent interest 
in exploring such models (e.g. 0, 13, IE H 0, Q)- To ex- 
plain the glitches as other than statistical flukes, these 
models require order unity variations in the curvature 
power spectrum across about an e-fold in wavenumber. 

Such cases are typically handled by numerically solv- 
ing the field equation on a case- by-case basis (e.g. [§]). 
For model-independent constraints and model building 
purposes it is desirable to have a simple but accurate pre- 
scription that relates features in the inflaton potential to 
features in the power spectrum (cf. [lol [Til fl3t ) . 

The generalized slow roll (GSR) approximation was in- 
troduced by Stewart [l4[ to overcome some of the prob- 
lems of the ordinary slow roll approximation for poten- 
tials with small but sharp features. In this approxima- 
tion, the ordinary slow roll parameters are taken to be 
small but not necessarily constant. In this paper we ex- 
amine and extend the GSR approach for the case of large 
features where the slow-roll parameters are also not nec- 
essarily small. 

In CTfl we review the GSR approximation and develop 
the variant for large power spectrum features. In the 
Appendix, we compare this variant to other GSR ap- 
proximations in the literature [HI, EU, EH, EH, E3- We 
show that our variant provides both the most accurate 
results and is the most simply related to the inflaton 
potential. In £)III1 we show how this technique can be 
used to develop alternate inflationary models to explain 
a given observed feature. We discuss these results in ^IVI 



II. GENERALIZED SLOW ROLL 

The GSR formalism was developed to calculate the 
curvature power spectrum for inflation models in which 
the usual slow roll parameters, defined in terms of time 
derivatives of the inflaton field cf> and the expansion rate 
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are small but t]h{= —Si) is not necessarily constant. In 
these models, the third slow-roll parameter 
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can be large for a small number of e- folds [HI, EH, EH]- 
Here and throughout we choose units where the reduced 
Planck mass (SttG)" 1 / 2 = 1. 

We study here the more extreme case where r\n is also 
allowed to become large for a fraction of an e-fold. These 
models lead to order unity deviations in the curvature 
power spectrum. As we shall see, different implementa- 
tions of the GSR approximation perform very differently 
for such models. 

An example of such a case is a step in the inflaton 



potential of the form V(<p) = m^ s (0)^/2, where the 
effective mass of the inflaton potential is given by @ 



off 



(0) 



1 + c tanh 



(3) 



This form for the potential has been shown to be a good 
description of large features in the temperature power 
spectrum at I ~ 20 — 40 tentatively seen in the WMAP 
data maximum likelihood (ML) parameters 

values for WMAP 5 are b = 14.668, c = 1.505 x 10~ 3 , 
d = 0.02705 and m = 7.126 x 10" 6 @. The potential 



14.5 



14.6 



14.7 



14.8 



0/M p 



0.00 



V [Gpc] 



FIG. 1: Upper panel: inflationary potential with a step from 
Eq. <j3j with parameters that maximize the WMAP5 likeli- 
hood (ML, black/solid) and an m 2 <p 2 potential that matches 
the WMAP5 normalization (smooth, red/dashed). Lower 
panel: conformal time to the end of inflation as a function 
of the value of the field. 



for this choice of parameters is shown in Fig. [T] (upper 
panel) . For comparison we also show the best fit smooth 
model (c = 0) with m = 7.12 x 1CT 6 . Since it will be 
convenient to express results in terms of physical scale 
instead of field value, we also show in the lower panel the 
relationship to the conformal time to the end of inflation 
7] = J t * ond dt'/a. Note that rj is defined to be positive 
during inflation. The two models have comparable power 
at wavenumbers k ~ ?y _1 ~ 0.02 Mpc -1 . 

The slow-roll parameters for these models as a function 
of •>] arc shown in Fig. [51 Notice that e# remains small 
in the ML model though its value changes fractionally by 
order unity. On the other hand, tjh is of order unity and 
82 is greater than unity in amplitude in this model around 
77 ~ 1 Gpc when the inflaton rolls across the feature. 



FIG. 2: Slow-roll parameters en, f]H and 52 for the two models 
of Fig. [l] ML step model (black/solid) and smooth model 
(red/dashed). 



Following [3], we begin by transforming the field equa- 
tion into dimensionless variables y = \/2kUk, x = krj 



dx 2 \ x 2 1 ' x 2 



where 
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Primes here and throughout arc derivatives with respect 
to In 77. 

The functions / and g carry information about devia- 
tions from perfect slow roll e# = 0, t]h = and 5 2 = 0. 
Specifically, without assuming that these three parame- 
ters are small or slowly varying 



A. Exact Relations 

It is useful to begin by examining the exact equations 
and solutions. The exact equation of motion of each k- 
mode of the inflaton field is given by 



d 2 u k 
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where 
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The field amplitude is related to the curvature power 
spectrum by 



A 2 



k 3 

(k) = —r lim 
27T 2 ki^a 



Uk 
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f 

— = -aHr)(e H - r) H ) + (1 - aHrj) , 
j = 3^ + 2[(aH V ) 2 -l} 

+(aHr]) 2 [26H - + 2e 2 H - Arj H e H + S 2 ] , 



(9) 



and the dynamics of the slow-roll parameters themselves 
are given by 



de H 
din a 
dr\n 
din a 



= 2e H (e H - T] H ) , 
= chVh + Vh - $2 ■ 



(10) 
(11) 



Moreover, these quantities are related to the inflaton po- 
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FIG. 3: Source functions for the deviations from slow roll 
used in the GSR approximations: 2g/3, 2gv/3 and G' (see 
§11 C[) for the maximum likelihood model. To good approxi- 
mation g — gv which directly relates the source function to 
features in the inflaton potential. Likewise G' ~ 2gv/3 and 
is most simply related to the curvature power spectrum for 
large deviations. 



tential by 
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1 - e H /3 



(12) 



which in the limit of small and nearly constant tjh, £h 
return the ordinary slow roll relations. 

In general, there is no way to directly express the 
source function g in terms of the potential without ap- 
proximation. Here we want to consider a situation where 
the feature in the potential is not large enough to inter- 
rupt inflation and hence tu <C 1, but is sufficiently large 
to make t]h of order unity for less than an e-fold. By 
virtue of Eq. (fTTj) . \S 2 \ ^> 1 during this time. This differs 
from other treatments which assume \ vh\ <C 1 and by 
virtue of Eq. ([TTJ]) a nearly constant e# [l4[ . 

Even under these generalized assumptions there are 
some terms in tjh and 82 that can be neglected. For 
example, even if t\h is not small, it suffices to take 



aHi-j - 1 = e H + e H 0(i] H ) 



(13) 



This expression preserves the ordinary slow roll relations 
when \t]h\ <C 1. When tjh is not small, this quantity 
remains of order en and so is negligible compared with 
bare rjn and 62 terms. Hence this approximation suffices 
everywhere. 

Following this logic, we obtain 



g = g v + e H 0(e H ,'q H ,52) , 
where gv is directly related to the potential 
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e H 0{e H ,r)H,6 2 ). (15) 



As shown in Fig. [3l this relationship between the source 
function g and features in the potential V holds even for 
the ML step potential. Thus, if we can express the func- 
tional relationship between g and the curvature power 
spectrum that is valid for large g we can use features in 
the power spectrum to directly constrain features in the 
inflaton potential. 

To determine this relation note that in the x and y 
variables the curvature is 1Z — xy/f, and its power spec- 
trum is A^(fc) = lim x< i \R\ 2 . The LHS of Eq. © is 
simply the equation for scale invariant perfect slow roll 
and is solved by 



(16) 



y (x)= 1 + - e", 



and its complex conjugate y$(x). An exact, albeit formal 
solution to the field equation can be constructed with the 
Green function technique [lj| 



y(x) = Vo(x) 



du 



g (In u)y {u)lm[y q (u)y o (x)} . (17) 



The solution is only formal since y appears on both the 
left and right hand side of the equation. The correspond- 
ing formal solution for the curvature power spectrum can 
be made more explicit by parameterizing the source y(u) 
as 



y(u) = F R (u)Rc[y (u)} + iFi(u)lm[y (u)j 
so that 



lim (xy) 

K<1 



i- \ I —^Fi{u)g(hxu) 



(18) 



(19) 
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du 
u 

du 
u 

c du 
u 



— W(u)Fi(u)g(\nu) 
X{u)F R (u)g{\nu) 
W(u)F K (u)g{lnu) + 0(x 2 ), 



where 



W{u) 



X(u) 



Im[y (u)]Rc[y (u)] 

u 

3sin(2u) 3cos2w 3sin(2w) 
2u 3 u 2 2u 

3 

-Re[y (u)]Re[y (u)] 
u 

3cos(2u) 3sin(2u) 

2m 3 u 2 
3 cos(2u) 



2 a 



Note that \\m u ^ W(u) = 1 and \im u ^ X(u) 
and we have utilized the fact that 



(20) 
u 3 /3 



lm.[y (u)]lm[y (u)] = 1 



1 



-X{u) (21) 
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and so Eq. (jT9|) becomes 



1.2 



FIG. 4: Ratio of field solution y to the scale invariant ap- 
proximation yo- Upper panel: real part Fr for a smooth case 
(red/dashed line), and for the maximum likelihood model 
(black/solid line), both at k — 10 -4 Mpc -1 . Lower panel: 
imaginary part Fi for the same models. 



goes to 1/u 2 in the limit u —* 0. 

Finally, the curvature power spectrum becomes 

Jlm(y)] 2 + [Rc(y)] 2 



Al(k) = lim x 2 

X<1 



f' 



(22) 



with y given by Eq. (| 1 9|) . 



lim 



T^gsrsI — j 



it 
3/ 



(24) 

where note that we have dropped the Re(xy) contribution 
since it adds in quadrature to the power spectrum and 
hence is second order in g. We call this the "GSRS" 
approximation for the curvature power spectrum A 2 ^ = 
lim^i |7?-gsrs| 2 given its validity for small fluctuations 
in the field solution from y — > yo- 

The choice of x is somewhat problematic (T^j. From 
Fig. [4] we see that taking x too small will cause spurious 
effects since Fi increases as x decreases. On the other 
hand, x cannot be chosen to be too large for the ML 
model since it will cause some k modes to have their cur- 
vature calculated when the inflaton is crossing the fea- 
ture. Moreover, if x is set to be some fixed conformal 
time during inflation ry m i n , then it will vary with k. 

We illustrate these problems in Fig. [51 For ?7 m i n = 10" 1 
Mpc (upper panel), GSRS undcrprcdicts power at low 
k for the smooth model and overpredicts it for the ML 
model. Agreement for the smooth model is improved 



by choosing 



10- 



i.e. nearer to horizon crossing 



B. GSR for Small Deviations 



(cf. Appendix for variants that take x w 1). On the 
other hand, the agreement for the ML model becomes 
worse and has a spurious feature at k ~ 10~ 5 Mpc -1 
where the inflaton is crossing the feature at x = 10~ 2 . 
In the next section, we shall examine the origin of the 
deviations from the exact solution and how a variant of 
the GSR approximation can fix most of them. 



The fundamental assumption in GSR is that one re- 
covers a good solution by setting F\{u) = Fr(u) = 1 in 
the formal solution for the field fluctuations in Eq. (fT9|) . 
Equivalcntly, y(u) — ► ya(u) in the source term on the 
RHS of Eq. ||7J). Note that this does not necessarily re- 
quire that g itself is everywhere much less than unity. 
For example, modes that encounter a strong variation in 
g while deep inside the horizon do not retain any im- 
print of the variation and hence the GSR approximation 
correctly describes the curvature they induce. 

In Fig. |U we show an example of F\ and Fr for a 
mode with k = 10~ 4 Mpc -1 for both the ML and smooth 
models. For the ML model, this mode is larger than 
the horizon when the inflaton crosses the feature. Note 
that even in the smooth model, the two functions deviate 
substantially from unity at x <C 1. In fact, they continue 
to increase indefinitely after horizon crossing and Fr oc 
x~ 3 diverges to compensate for |Re(j/o)| °^ x 2 . For the 
ML model, even Fi deviates strongly from unity during 
the crossing of the feature at x ~ 0.1. 

The impact that these deviations have on the curvature 
spectrum can be better understood by reexprcssing the 
various contributions in a more compact form. First note 
that 

,. a; 3 f°° du , „ , 1 ff'\ 

lim - / — u~ 3 g (Inu = "« T ' 23 

z«i 3 J x u 3 V / J 



C. GSR for Large Deviations 

When considering large deviations from scale invari- 
ance, either due to sharp features in the potential or 
due to extending the calculation for many e-folds after 
horizon crossing, the first qualitative problem with the 
GSRS approximation of Eq. (|2"4")) is that it represents a 
linearized expansion for a correction that is not necessar- 
ily small. When the correction becomes large, 7?-gsrs can 
pass through zero leaving nodes in the spectrum. While 
this is not strictly a problem for the ML model, it is bet- 
ter to have a more robust implementation of GSR for 
likelihood searches over the parameter space. 

We can finesse this problem by replacing the linearized 
expansion 1 + x by e x and write the power spectrum in 
the form 

2 r°° dn 

\nA 2 n (k) = G(\nn nlin ) + - / -^(fc^^ln^) , (25) 
where 

G(ln»/)=ln^+H^. (26) 

This procedure returns the correct result at first order 
since g and /'// are both first order in the slow-roll pa- 
rameters (see Eq. ©). We shall sec below that it can be 
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FIG. 5: GSRS approximation to the curvature power spec- 
trum (dashed lines) compared to the exact solution (solid 
lines) for a choice of ?7 m in = 1CP 1 Mpc (upper panel) and 
?7mm = 10~ 2 /fc Mpc (lower panel). The ML model is shown 
in blue and the smooth model in red for GSRS. 



further modified to match the fully non-linear result for 
superhorizon modes. 

The more fundamental problem with GSRS is the devi- 
ation of the true solution y from the scale invariant solu- 
tion yo when the mode is outside the horizon (see Fig.QI. 
The origin of this problem is that the exact solution re- 
quires the curvature 1Z = xy/f to be constant outside the 
horizon, independently of how strongly / evolves. Thus, 
if / is allowed to vary significantly, either due to the 
large number of e-folds that have intervened since hori- 
zon crossing or due to a feature in the potential, then y 
must follow suit and deviate from yo breaking the GSRS 
approximation. 

Fig. [6] (upper panel) illustrates this problem. Even for 
the smooth model, the curvature is increasingly under- 
estimated as x — > . With the ML model, the crossing 
of the feature induces an error of the opposite sign. For 
x ~ 1CP 5 these problems fortuitously cancel but not for 
any fundamental or model independent reason. 

Given this problem, GSRS actually works better than 
one might naively expect. For example at k = 10~ 4 
Mpc" 1 , even though F l ~ 1.28 at x = 1CT 5 , the GSRS 
approximation gives a ~ 2.5% difference in the curvature 
and a ~ 5% difference in the power spectrum with the 
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FIG. 6: Curvature evolution after horizon crossing in the 
GSRS (upper panel) and GSRL (lower panel) approxima- 
tions, both normalized to the exact solution. The ML model 
(black/solid line) and smooth model (red/dashed line) are 
both shown at k = 1CT 4 Mpc -1 . 



exact solution for the smooth model instead of the 28% 
and the 64% differences one might guess. The main con- 
tribution to the GSRS correction from scale invariance is 
given by the integral term in Eq. (f24|) , which is ~ 0.25 
for the smooth case. Given that Fi is a linear function 
in In 77 and g is slowly varying, we can approximate en- 
hancement due to Fi of the integral term by its average 
interval (~ 1.14). With this rough estimate we obtain an 
approximately (1 + 0.25) 2 /(l + 0.25 x 1.14) 2 - 5 - 6% 
error in power in agreement with the power spectrum 
result in Fig. 

Furthermore although Fr diverges as x~ 3 in Fig. 31 
the contribution to the power spectrum of the real part 
of y remains small. Its absence in the GSRS approxima- 
tion produces a negligible effect for modes that are larger 
than the horizon when the inflaton crosses the feature. 
The integrands for the real contribution contain either 
the function X 7 which peaks at horizon crossing x ~ 1, 
or x W(u) which is likewise suppressed at x <C 1. The 
correction adds in quadrature to the imaginary part and 
so it is intrinsically a second order correction (see £ 111 El) . 
For k = 10~ 4 Mpc -1 its contribution to the power spec- 
trum is 0.08% of the power spectrum in the ML model. 
The fact that integrals over the deviation of y from yo 
can remain small even when neither g nor the maximum 
ofy — yo is small is crucial to explaining why the GSR ap- 
proximation works so well and why we can extend GSRS 
with small, controlled corrections. 

Nonetheless these problems with GSRS are significant 
and exacerbated by the presence of sharp features in the 
potential. The fundamental problem with GSRS is that 
its results depend on an arbitrarily chosen value of x <C 1 , 
i.e. 1Z is not strictly constant in this regime. Phrased in 
terms of Eq. (|2"6"]l the problem is that g is not directly 



6 



related to G but rather 
2 



3- 9 = G ' + 3 



where 



G' = 



dG = 2 f" f f 2 



(27) 



(28) 



In GSRS, replacing g with 3G'/2 amounts to a second 
order change in the source function. In fact even for 
the ML step function this change is a small fractional 
change of the source everywhere in In 77: it is small as the 
inflaton rolls past the feature since |/"//| 3> (f'/f) 2 and 
it is small before and after this time since ^ 1- I n 

terms of the slow-roll parameters, this replacement is a 
good approximation if r/j^ ~ 0(1) only where l^l 3> 1 
and g w S 2 (see Eqs. (14]) and (fTS"]) ). 



(29) 



ni ^ 2 2 
3 3 



e H 0{e H ,r]H,S2 



Moreover G" « 2gy/3 and remains directly relatable to 
the inflaton potential through Eq. (|15p . For comparison 
we show all three versions of the GSR source function in 
Fig.H 

Nonetheless, the replacement can have a substantial 
effect on the curvature once the source is integrated over 
ln?7 because the difference is a positive definite term in 
the integral. Moreover, this cumulative effect is exactly 
what is needed to recover the required superhorizon be- 
havior. Replacing 2^/3 — s- G' in the power spectrum 
expression, we obtain [13| 



lnA^(fc) = G(ln77„ 



^-W(kri)G'(]nri), (30) 

Vrpln V 



which we call the GSRL approximation. The field so- 
lution corresponding to this approximation, valid for 
x <C 1, is given by 



lim \xy\ = exp 

ai<gl 



1 /' 1 

3j + 2 



fill 

— W(u)G'(\nu) 
u 



(31) 



Now any variation in / while the mode is outside the 
horizon and W(krj) pa 1 integrates away and gives the 
same result as if In ?7 m i n were set to be right after horizon 
crossing for the mode in question. This can be seen more 
clearly by integrating Eq. (f30|) by parts [l]| 



lnA^(fc) 



ch) 



W'(krj)G(lnT]). (32) 



Since — J °° d\nxW'(x) = 1 and lim^o W (x) — > 0, the 
curvature spectrum does not depend on the evolution of 
/ outside the horizon. Moreover, the integral gets its 
contribution near x ~ 1 so for smooth functions G(hir]) 
we recover the slow roll expectation that 



InA^(ft) « G(lnr?) 



(33) 
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FIG. 7: GSRL approximation to the curvature power spec- 
trum. Upper panel: approximation compared with the ex- 
act solution (solid lines) for the maximum likelihood model. 
Lower panel: fractional error between the approximation and 
the exact solution. 



If the slow-roll parameters are all small then the leading 
order term in Eq. (|33|) returns the familiar expression 



for the curvature spectrum w f~ 2 w H 2 /8it 2 eh at 
krj « 1. Choe et al. [l5[ showed that Eq. (|32|) is correct 
up to second order in g for krj <^ 1. Here we show that it 
is correct for arbitrary variations in / and g outside the 
horizon. 

The superhorizon curvature evolution for k = 10~ 4 
Mpc -1 corresponding to the GSRL approximation is 
shown in Fig. [6] (lower panel). In the x <C 1 domain of 
applicability of Eq. pip , the curvature is now appropri- 
ately constant for both the ML and smooth models. The 
net result is that the curvature power spectrum shown in 
Fig. [7] is now a good match to the exact solution for low 
k. 



D. Power Spectrum Features 

We now turn to issues related to the response of the 
field and curvature for k modes that are inside the hori- 
zon when the inflaton rolls across the feature. Fig. [7| 
shows that the GSRL approximation works remarkably 
well for the ML model despite the fact that the power 
spectrum changes by order unity there. The main prob- 
lem is a ~ 10 — 20% deficit of power for a small range in 
k near the sharp rise between the trough and the peak. 

In Fig. [8j we show the deviation of the exact solution 
y from the scale invariant yo that is at the heart of the 
GSR approximation. The three modes shown, fc^ip = 

1.8 X 10~ 3 MpC~\ fc nodc = 2.5 X 10~ 3 MpC~\ fcbump = 

3.2 x 10 -3 Mpc -1 , correspond to the first dip, node and 
bump in the power spectrum of the ML model. 
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FIG. 8: Fractional difference between \yo\ and \y\ for the ML 
model at k values at the dip, node and bump of the feature 
in the power spectrum (see text). 



The first thing to note is that for higher k, the inflaton 
crosses the feature at increasing x where the deviations 
of y from yo actually decrease. Hence the fundamen- 
tal validity of the GSR approximation actually improves 
for subhorizon modes. Combined with the GSRL ap- 
proximation that enforces the correct result at x <C 1, 
this makes the approximation well behaved nearly every- 
where. 

The small deviations from GSRL appear for modes 
that cross the horizon right around the time that the 
inflaton crosses the feature. It is important to note that 
the step potential actually provides two temporal fea- 
tures in g or G' displayed in Fig. [3] Each mode first 
crosses a positive feature at high r\ and x and then goes 
through a nearly equal and opposite negative feature. 
The end result for the field amplitude or curvature is an 
interference pattern of contributions from both temporal 
features. For example, the peak in power is due to the 
constructive interference between a positive response to 
the positive feature and a negative response to the neg- 
ative feature. This suggests that one problem with the 
GSRL approach is that it docs not account for the de- 
viation of the field y from yo that accumulates through 
passing the positive temporal feature when considering 
how the field goes through the negative feature. This is 
intrinsically a non-linear effect. 

The final thing to note is that, since g and G' are of 
order unity as these modes exit the horizon, the real part 
of the field solution is not negligible. Moreover, it con- 
tributes a positive definite piece to the power spectrum. 
In Fig. [9l we show the result of dropping the real part 
from the exact solution. Note that the fractional error 
induced by dropping the real part looks much like the 
GSRL error in Fig. [7] but with ~ 1/2 the amplitude. 
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FIG. 9: Contribution of the real part of the y field to the 
curvature power spectrum. Upper panel: spectrum with and 
without the real part. Lower panel: fractional error between 
the two solutions. 



E. Iterative GSR Correction 

The good agreement between GSRL and the exact so- 
lution even in the presence of large deviations in the cur- 
vature spectrum suggests that a small higher order cor- 
rection may further improve the accuracy. Moreover, the 
analysis in the previous section implies that there are two 
sources of error: the omission of the field response from 
inside the horizon x > 1 when computing the response 
of the field to features at horizon crossing x ~ 1 and the 
dropping of the real part of the field solution. 

Both of these contributions come in at second order in 
the GSR approximation. All first order GSR variants in- 
volve the replacement of the true field solution y with the 
scale invariant solution yo in Eq. (|7|). This replacement 
can be iterated with successively better approximations 
to y. We begin with the GSRS approximation of replac- 
ing y — > yo to obtain the first order solution y\ . We then 
replace y — > y\ in the source to obtain a second order 
solution y2, etc. 

We show the fractional error between the iterative so- 
lutions and the exact solution for k = fcnode i n Fig. 1 1 01 
where the error in GSRL is roughly maximized. As in 
the first order GSRS approach, the accuracy depends on 
the arbitrary choice of x = krj m i n when the curvature is 
computed. The number of iterations required for a given 
accuracy increases with decreasing x. We show the cur- 
vature spectrum in Fig. [IT] for the same two choice of 
Vmm = 10 _1 Mpc (upper panel) and ry m i n = 10~ 2 //c Mpc 
(lower panel) as in Fig. [5] Note that in both cases, the 
result has converged at the 0.5% percent level or better 
to the exact solution within three iterations. 

Unfortunately the iterative GSRS approach is not of 
practical use in that each iteration requires essentially 
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FIG. 10: Fractional difference between the exact (y) and nth 
order iterative solutions (y n ) for the ML model at k — fc no d c 
where the errors in the GSRL approximation are maximized. 
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FIG. 11: Curvature power spectrum in the GSRS approxima- 
tion for r)xnm = 10 -1 Mpc (upper panel) and n m i n = 10~ 2 /fc 
Mpc (lower panel) when y — > y n in the GSRS source com- 
pared to the exact solution. 



the same effort as a single solution of the exact approach. 
On the other hand, rapid convergence in the iterative 
GSRS approach suggests that a nonlinear correction to 
GSRL based on a second order expansion might suffice. 
A second order GSRL approach differs conceptually from 
the iterative GSRS approach in that it is formally an 



expansion in g where in our case \g\ <C 1 is not satisfied. 
The iterative GSRS approach is exact in g but expands 
in y — yo- What makes a second order GSRL approach 
feasible is that the critical elements involve time integrals 
over g which can be small even if g is not everywhere 
small. 

Our strategy for devising a non-linear correction to 
GSRL is to choose a form that reproduces GSRL at first 
order in g, is exact at second order in g, is simple to relate 
to the inflaton potential, and finally is well controlled at 
large values of g. The second order in g expressions for 
the curvature are explicitly given in [loT ] and come about 
by both iterating the integral solution in Eq. (fT5|) and 
dropping higher order terms. Our criteria are satisfied 
by 



1 



1 



A4 = A^| GSRL [l + -/ 1 2 (fc) + -/ 2 (fc)] 



where 



(34) 



h{k) 
Hk) 



1 


ran 


71 


Jo V 


-4 


c * 




'o u 



with 



F{u) 



dvf 

V 2 f 



(35) 



(36) 



We call this the GSRL2 approximation. In the Appendix 
we discuss alternate forms [r| . 

In the GSRL2 approach, I\ corrections come half from 
the first order calculation of the real part of the field and 
half from iterating the imaginary part to second order. In 
Fig. [12] we show 1\ and I2 for the ML model. Note that 
If dominates the correction to the net power as it always 
enhances power, while I2 is both smaller and oscillates in 
its correction. Furthermore, both <C 1 and I/2I <C 1 
for the ML model which justifies a second order approach 
to these corrections. The GSRL2 correction can be taken 
to be {1 + If + 1 2} in this limit. 

We show in Fig. [Tj|] how the GSRL2 corrections re- 
duce the power spectrum errors of GSRL in Fig. [7] for 
the ML model. For the full GSRL2 expression the power 
spectrum errors are reduced from the 10 — 20% level to 
the < 4% level. We show that the GSRL2 approxima- 
tion remains remarkably accurate for substantially larger 
features in the Appendix. 

Moreover, the errors are oscillatory and their observ- 
able consequence in the CMB is further reduced by pro- 
jection. The temperature and polarization power spectra 
are shown in Fig. [14] and [15] and the errors are < 0.5% 
and < 2% for the respective spectra. 

Given the intrinsic smallncss of I2 and its oscillatory 
nature, the most important correction comes from the 
positive definite I\ piece. Note that it is a single integral 
over the same G function as in the linear case. Thus, 
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FIG. 12: Second order GSRL2 power spectrum correction 
functions if and Ii for the ML model. 
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I\ corrections simply generalize the GSRL mapping be- 
tween G and curvature in a manner that is equally simple 
to calculate. I2 on the other hand is more complicated 
and involves a non-trivial double integral with a different 
dependence on the inflaton potential. 

We also show in Figs. IT3IT51 the results for the GSRL2 
expression with I2 omitted. While the curvature power 
spectrum errors increase slightly to ~ 5%, the temper- 
ature power spectrum errors at < 2% are still well be- 
low the ~ 20% cosmic variance errors per I at I ~ 30. 
They are even sufficient for the cosmic variance limit of 
coherent deviations across the full range of the feature 
(20 < I < 40) 20%/v / 20 ~ 4 - 5% in the ML case. The 
polarization spectrum has slightly larger errors due to 
the reduction of projection effects but still satisfies these 
cosmic variance based criteria. 



III. APPLICATIONS 

In the previous section, we have shown that a par- 
ticular variant of the GSR approximation which we call 
GSRL2 provides a non-linear mapping between G' and 
the curvature power spectrum. G' quantifies the devi- 
ations from slow roll in the background and moreover 
is to good approximation directly related to the inflaton 
potential. These relations remain true even when the 
slow-roll parameter tjh is not small compared to unity 
for a fraction of an e-fold. 

This relationship is useful for considering inflation- 
model independent constraints on the inflaton potential. 
It is likewise useful for inverse or model building ap- 
proaches of finding inflaton potential classes that might 
fit some observed feature in the data. We intend to fur- 
ther explore these applications in a future work. 

Here as a simple example let us consider a poten- 
tial that differs qualitatively from the step potential but 
shares similar observable properties through G': V(<j)) = 



FIG. 13: GSRL2 approximation to the curvature power spec- 
trum. Upper panel: approximation of Eq. (|34p (red/dashed 
line) compared to the exact solution (black/solid line). We 
also show the GSRL2 approximation omitting the I2 term 
(blue/dashed-dotted line). Lower panel: fractional error be- 
tween these GSRL2 approximations and the exact solution. 
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FIG. 14: GSRL2 approximation to the CMB temperature 
power spectrum. Upper panel: approximation (red/dashed 
line) compared to the exact solution (black/solid line). We 
also show the GSRL2 approximation omitting the I2 term 
(blue/dashed-dotted line). Lower panel: fractional error be- 
tween the GSRL2 approximations and the exact solution. 



a transient perturbation instead of a step 



'off 



1 + Ae -(<t>-bf/(2* 2 )( (j) _ b) 



(37) 



'err 1 ' 



/2 where the effective mass of the inflaton now has 



In Fig. [16] we show the potential for the choice of pa- 
rameters b = 14.655, A = 0.0285, a = 0.025, and 
m = 7.126 x 10 -6 (upper panel) and we also show G' 
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FIG. 15: GSRL2 approximation to the CMB _E-mode polar- 
ization power spectrum. The same as in Fig. 1141 
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FIG. 16: Alternate inflationary model with a perturbation in 
the mass. Upper panel: comparison of potential in Eq. (|37p 
(black/solid line) and the smooth potential (red/dashed line). 
Lower panel: source function of the deviation from slow roll 
G' for the same models. 



in the lower panel. For comparison we show the smooth 
case A = 0. Comparison with Figs. Q] and [5] shows that 
this potential, which has a bump and a dip instead of 
a step, produces a similar main feature in G but has 
additional lower amplitude secondary features. 

In Fig. [T7]we compare the GSRL2 approximation with 
and without the double integral I2 term compared to the 
exact solution. Notice that GSRL2 performs equally well 
for this very different sharp potential feature. Further- 
more, similarity in G with the step potential carries over 
to similarity in the curvature power spectrum. 



FIG. 17: GSRL2 approximation to the alternate model of 
Fig. 1161 Upper panel: approximation (red/dashed line) com- 
pared to the exact solution (black/solid line) for an effective 
mass given by Eq. (|3T[) . We also show the GSRL2 approx- 
imation with I2 omitted (blue/dashed-dotted line). Lower 
panel: fractional error between GSRL2 approximations and 
the exact solution. 



IV. DISCUSSION 

We have shown that a variant of the generalized slow 
roll (GSR) approach remains percent level accurate at 
predicting order unity deviations in the observable CMB 
temperature and polarization power spectra from sharp 
potential features. Unlike other variants which explic- 
itly require \r\n\ «C 1, and hence nearly constant e#, our 
approach allows i]h to be order unity, as long as it re- 
mains so for less than an e-fold, and hence ch to vary 
significantly. We have tested our GSR variant against 
a step function model that has been proposed to ex- 
plain features in the CMB temperature power spectrum 
at i ~ 20 - 40. 

Our analysis also shows that to good approximation a 
single function, G'(ln?7), controls the observable features 
in the curvature power spectrum even in the presence of 
large features. We have explicitly checked this relation- 
ship and the robustness of our approximation by con- 
structing two different inflationary models with similar 
G . 

Therefore observational constraints from the CMB can 
be mapped directly to constraints on this function in- 
dependently of the model for inflation. Moreover, this 
function is also simply related to the slope and cur- 
vature of the inflaton potential in the same way that 
scalar tilt is related to the potential in ordinary slow roll 
G « 3(V,p/V) 2 - 2{V^/V). These model independent 
constraints can then be simply interpreted in terms of 
the inflation potential. We intend to explore these appli- 
cations in a future work. 
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FIG. 18: Curvature power spectrum for the ML and 3ML 
models. 
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APPENDIX A: OTHER GSR VARIANTS 

In this Appendix, we compare various alternate forms 
discussed in the literature for the curvature power spec- 
trum under the GSR approximation. We test these ap- 
proximations against the GSRL and GSRL2 approxima- 
tions of the main text for the ML model and a more 
extreme case with c = 3cml = 0.004515 (with other pa- 
rameters fixed) denoted 3ML (see Fig. [T5)) . We begin 
by considering variants that are linear in the GSR ap- 
proximation and then proceed to second order iterative 
approaches. 

The first variant is the original linearized form of GSRS 
given in [3 ("S02") 



1 



1 



2f 
3 / 



du 



W(u)g(\Yiu) 



(Al) 

Like GSRS, this approximation depends on an arbitrary 
choice of x but its impact is exacerbated by the lin- 
earization of the correction here. In Fig.QlJlwe show the 
fractional error in this approximation for rj m i n = 10 _1 
Mpc. Note that because of the linearization, the curva- 
ture power spectrum can reach the unphysical negative 
regime (shaded region). 

A second variant further exploits the relationship be- 
tween the GSR source functions /, /'// and g and the po- 
tential through the slow-roll parameters (see Eq. ([9])). By 
further assuming that \t]h\ <C 1, terms involving V^/V 
can be taken to be constant and evaluated instead at hori- 
zon crossing k = aH (sec Eq. (HHJ))- Finally by rewriting 



FIG. 19: Fractional error in the curvature power spectrum 
for first order GSR variants for the ML model (lower) and 
the 3ML model (upper). 



the change in /'// as the integral of (/'//)', one obtains 
[3 ("DS02") 



-2 



,<t> 
du 



(A2) 



where a w 0.73 and with r/ ss 1/aH, u = k/aH. Here 



W e (u*,u) = W(u)-9(u*-u) 



(A3) 



with the step function 9{x) = for x < and 9(x) = 1 
for x > 0. Note that lim„_»o W$(l, u) = and hence the 
function has weight only near horizon crossing at u 1. 

For cases like the ML and 3ML models where r/n is nei- 
ther small nor smoothly varying, these DS02 assumptions 
have both positive and negative consequences. They 
largely solve the problem for superhorizon modes dis- 
cussed in ^11 CI by extrapolating the evaluation of the 
potential terms from kr\ <C 1 to kr\ ~ 1. On the other 
hand, a large j]h means that en evolves significantly. 
Artifacts of this evolution appear through the prefactor 
(V/V i( p) 2 cx 1/ch in Eq. (|A2|) most notably in the form 
of a spurious feature at k ~ 10 ~ 3 Mpc -1 in Fig. [T91 
Finally, like S02, DS02 does not guarantee a positive def- 
inite power spectrum. 
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FIG. 21: Temperature power spectrum for c = 8cml (and the 
other parameters fixed). Note that even in this extreme, ob- 
servationally unviable, case the temperature power spectrum 
has < 22% errors everywhere for GSRL2 whereas the linear 
GSRL and CGS04c approximations substantially under and 
over predict power respectively. 



x* to obtain ("CGS04a") 



FIG. 20: Fractional error in the curvature power spectrum 
for second order GSR variants (see text) for the ML model 
(lower panel) and the 3ML model (upper panel). Note that 
the error in CGS04a has been divided by a factor of 10 for 
plotting purposes. 



A third variant is to replace G' with 2gv/3 in Ec 
([3TJ]) so that the source directly reflects the potential [R 
("KDHS05") 



lnA^(A;) = G(lnr] min ) + 



W{krj)g v • (A4) 



A 2 



U 1 

+ ( 3a 
1 



2a 



J * 



5^ 
~L2~ 



—at 
3 * 



12/ /* 



3a 



2 

12 L "^ 
7 



8a* + — 7r 2 a* 



f* 
2C(3) 



-.4 



f * J * 



(A5) 



As we have seen in SjTTC] this approximation is actually 
fairly good locally in In 77 and hence locally in k around 
the feature. However the omission of rj^j terms causes a 
net error in the spectrum for k modes that cross out of the 
horizon before the inflaton reaches the feature. Hence like 
the GSRS approximation, KDHS05 ovcrpredicts power 
at low k for the ML and 3ML models. Fig. [TO] shows a 
choice with ^ min = 10 _1 Mpc. 

We consider next second order GSR variants. The first 
variant [l5| begins with a second order approach as in 
GSRL2 but then further assumes that functions such as 
/'// carL be approximated by a Taylor expansion around 



where A = -4aJ + 16a* - 5/37r 2 a* - 8 + 6C(3), ( is 
the Ricmann zeta function, and a* = a — ln(a:*). This 
approach is essentially a standard slow roll approxima- 
tion carried through to third order with the help of an 
exact solution for power law inflation. For the ML and 
3ML models, applying this approximation leads to qual- 
itatively incorrect results as one might expect. We show 
this variant in Fig. [20] with x* = 1. 

A second variant attempts to retain both the generality 
of GSR and the evaluation of central terms at horizon 
crossing by implicitly modifying terms of order (/'//) 3 
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and higher when compared with GSRL2 ( "CGS04b" ) related to GSRL2 as integration by parts shows 



lnA^ 



In 



2/: , 1 (fl 



ft) 3/* 9 V/. 
— Wg\u*,u)g{u) 



(A6) 



du 



X{u)g(u) 



-- / —X(u)g(u) / —g(v) 

•J JO u Ju v 

du f°° dv 

— X g (u*,u)g(u) I -rg{v) 



where Wg was given in Eq. (|A3|) and 



Xg(u^,u) = X(u) —Q(u* — u) 



(A7) 



Here, the subscript * denotes evaluation near horizon 
crossing. In Fig. [20] we show the result with ti* = 1. 
Notably it performs worse than the first order GSRL ap- 
proximation for the 3ML model. 

Finally, the last variant considered takes (l5j 
("CGS04c") 



Jo w 



(A8) 



+2 

-4 



(A9) 



The main difference is that the second order corrections 
are exponentiated. This causes a noticeable overcorrec- 
tion for the 3ML model when compared with GSRL2. In 
Fig. [201 we compare the three variants mentioned above. 

Furthermore, in spite of the 20 — 40% errors in the 
curvature power spectrum in the 3ML model for GSRL2, 
the CMB temperature power spectrum has only 1 — 2% 
errors for £ > 20 and a maximum of < 5% errors at 
I < 20. As discussed in the text, this level of error is 
sufficient for even cosmic variance limited measurements 
at the I < 40 multipolcs of the feature. This reduction is 
due to the oscillatory nature of the curvature errors and 
projection effects in temperature. 

In fact for even larger deviations GSRL2 still performs 
surprisingly well for the temperature power spectrum. In 
Fig. [5T] we show the temperature power spectra for the 
GSRL2 approximation, and compare it with GSRL and 
CGS04c for a very extreme case with c = 8cml (and the 
other parameters fixed) . GSRL2 has a maximum of 22% 
error in the temperature power spectrum and predicts 
qualitatively correct features. Finally, the dominant cor- 
rection is from the term that is quadratic in 1\. The 
simplified GSRL2 form of 



1 K 



tcIgsrl 



(A10) 



works nearly as well. Thus, the curvature power spec- 
trum still depends only on G' to good approximation 
even in the most extreme case. 



where F(u) is given by Eq. (|3"r?)) . CGS04c is closely 
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